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1. INTRODUCTION 

In nonlinear dynamic systems, chaos synchronization is the first phenomenon which discovered by 
Fujisaka and Yamada in 1983, but did not receive great interest until 1990 when Pecora and Carrol developed 
this phenomenon between two identical chaotic systems with different initial condition [1-4]. Chaos 
synchronization has attracted considerable attention due to its important applications in physical systems [1], 
biological systems [5], Encryption [6] and secure communications [7], etc. After then, several attempts were 
made to create many types of synchronization phenomena such as Complete Synchronization (CS) [2, 4, 8], 
Anti-Synchronization (AS) [9, 10], Hybrid Synchronization (HS) [11], Generalized Synchronization (GS) 
[12], Projective Synchronization (PS) [13], Hybrid Projective Synchronization (HPS) [14] and Generalized 
Projective Synchronization (GPS) [15]. Amongst all types of synchronization schemes, PS and HPS attracted 
lots of attention because it can obtain faster communication in application to secure communication [13, 14]. 
Both of them are characterized that the two systems could be synchronized up to a constant diagonal matrix, 
but in the first feature, all diagonal elements of scaling matrix should be equal whereas these diagonal 
elements are different in the second feature. Obviously, choosing the constant matrix as unity will lead to CS. 
So, CS and AS are the special cases of PS and HS belong to the special case of hybrid projective 
synchronization. HPS gives more complexity to the controller and the message cannot be easily decoded by 
the intruder. 

In projective and HPS processes, various strategies have been introduced to stabilize dynamic error 
systems, including adaptive control [16], active control, nonlinear control [17-20] and linear feedback 
control [21-23]. Among many control strategies, the nonlinear control strategy has been continuously for 
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which great interest to many scientists, due to its effectiveness, reliability, and widely has been used as 
a single powerful strategy for synchronizing different class of the nonlinear dynamic systems [24, 25]. 
But, the control input design should be based on the functions of the controlled system according to the 
traditional nonlinear control. In order to simplify the control input, adaptive nonlinear control has been 
designed to facilitate the control input process. To ensure that the designed controller has a good control 
effect, the controller is designed for a non-linear control system based on the theory of stability Lyapunov 
with known and unknown parameters. Then, controllers designed to synchronize a hyperchaotic system were 
used. These findings may be important in understanding and controlling problems in modern society. 
Also the effectiveness and strength of controllers are verified by numerical simulation results. 


2. PROJECTIVE AND HYBRID PROJECTIVE SYNCHRONIZATION 

The PS and HPS are illustrated in this section. There are two nonlinear dynamical systems, the first 
is called drive system, whereas the second, is called response system, and the response system controls 
the drive system. The drive system and response system yields the (1) and (2), respectively [18] and [21]. 


x = Ax + f(x) (1) 
y=By+giy)+U (2) 


where A,B are a n Xn parameters matrices, x = [x,,X>,..,Xn]7 E R™ 14, y = [y1, Vo.) Val? E R™? are 
state vector, f(x) and g(y) are the nonlinear functions for system 1 and system 2, respectively. Also, 
U = [uy, Uz, ..., Un]” E R” is a control input vector. Whereas the error dynamical system is defined as 


e=y-—Sx (3) 


where e = [€),€3,...,€,]’ E R™? in general, S is n-order diagonal matrix i.e. S = diag(s1,S2,.--,Sy), S is 
called scaling matrix and $1, S2,...,Sp are called scaling factor. Our goal is to propose a suitable controller U 
to make the response system asymptotically approaches the drive system, and finally the synchronization 
phenomena will be achieved in the sense that the limit of the error dynamical system approaches zero i.e. 


lime ollell a limesolly > Sx|| =0 (4) 


The scaling matrix S play an important role to determine the phenomenon of synchronization, such 
as if S is constant matrix and 
— Sı = S2 = +++ = Sp, then this phenomenon is called PS 
— Sı É S2 #° £ Sp, then this phenomenon is called HPS 
— Ws; = 1, then this phenomenon is called CS 
— Ws; = —1, then this phenomenon is called AS 
— Ws; = +1, then this phenomenon is called HS 


3. APPLICATIONS 

In this section, we take 4-D non-linear dynamical system which discover by Zhang et al 
in 2017 [26], for example to show how to use the results obtained in this paper to analyse the synchronization 
class of hyperchaotic systems. The mathematical model is representing by the following: 


x=a(y-—x)- fw 
y=xz—aqy 
Z=b-xy-cz (5) 


w=ry—dw 


where x, y,z and w are state variable, xz and xy are two nonlinear terms and a, b,c,d, f,r and q are positive 
parameters. And the system (5) is hyperchaotic attractors due to has possess two positive Lyapunov 
exponents as LE, = 0.24, LE, = 0.23 at a=5, b = 20, c = 1, d = 0.1, f = 20.6, r = 0.1 andq = 1. 
The projections of hyperchaotic attractor of the above system are shown in Figure 1. 
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Figure 1. The attracters of the system (1) in: (a) y-w plane, (b) x-z-y space 


According to (1) and (2), system (5) can be represent as 


X1 -a a 0 —f]f*% 0 00 0 ofo 
ž%ļ|_|0 -q 0 0 |}%2 0 0 1 O Of}x,x3 
llo o -e ollel" lol tlo o -1 ollxx, (9 
X4 0 r 0 dll% 0 00 0 ollo 

A f (2) 
Yı —a a 0 -flM 0 00 0 ofo Uy 
Vo 0 -q 0 0|} 0 01 O Oftyiys Uz 
|= +|, |+ + 7 
V3 0 0 -c 0|» b 0 0 -1 Oflyiy2 U3 0) 
Va 0 r o0 -dally 0 00 o olo Ug 

B g) 


where U = [u, uz, Ug, u4] 7 is the controller to be designed 


4. PROJECTIVE SYNCHRONIZATION (PS) 
This phenomenon take place under the condition “that all diagonal elements s; of the constant 
scaling matrix S(t), possess the same value”. Herein, three cases are considered as 
- V Si = 3 
- Vs,=1 
- Vs,;=-1 


4.1. The Controllers at scaling factor V s; = 3 
According to (3), the error of PS é; € R* between the system (6) and the system (7) is depict by: 


êi = yi —3),x;, 1 = daig[1,1,...,1], i = 1,2,3,4 
and lead to 


é, = ale, — e1) — fey + u 
ê, = —qez + e3yı + 3(y1—x1)X3 + U2 


é3 = —Ce3 = e2y1 + 3 (xı — 1) X2 = 2b + U3 (8) 
é, = —de, + re, + uz 
Theorem 1. If design control U = [u,, uz, Uz, U4] 7 as: 
Uy = 0 
Uy = —ae, — re, + 3(x, — y1) x3 
_ E (9) 
uz = 2b + 304 — x1 )x2 
Uy = fey 
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Then the system (8) can be controlled i.e., PS between system (6) and system (7) is achieved. 
Proof: By inserting the controller (9) in the error system (8) we get: 


é, = ale, — e1) — fey 
ê, = —Gez — re, — Ae, + e3yı 


(10) 
€3 = —C&z — €2Y1 
È, => —de, + rez + fe, 
Construct the Lyapunov function as the following: 
tT =Ging(t ets 
V(e) =e’Pe, P = diag(>,5,5,5) (11) 


and derivative V (e) along time of (10) is: 


V(e) = €,(a(ez — e1) — fes) + €2(—Gez — re, — ae, + e3y1) 
+e3(—ce3 — e2y1) + e4(—de, + rez + fe) 


Above equation can reduce as: 
V(e) = —ae? — qe? — ce? — de? = —eTQ e, Q=diaglaqcd] 


the matrix Q is positive definite. So, V (e) is negative definite. The Lyapunov's direct method is satisfied. 
Therefore, the response system (7) is PS with the drive system (6) asymptotically, the proof is complete. 


4.2. The Controllers at scaling factor Y s; = 1 
For all scaling factor are equal 1, the error of PS between the drive system (6) and the response 


system (7) is given by: 


é, = ale, — e1) — fey + u 


Èz => —qe2 + €1X3 + €e3X1 +e, e3 + U2 (12) 
Èz = —Ce3 — €4Xz — 2X1 — €1€2 + U3 
é, = —de, + re, + u, 


Theorem 2. The systems (6) & (7) will be asymptotically stable, if the controller is designed as follows: 


Uy = —e3 X3 + €3X2 

Uz = —ae, — re, (13) 
uz =0 

u, = fey 


Proof: By substituting the controllers (13) in the system (12), we can obtain: 
é, = a(ez — e1) — fees — 2X3 + €3X2 
È> = —Gez + €1X3 + €3X1 + €1€3 — ae — re, (14) 
63 = —Ce3 — €4X2 — 2X1 — €1€2 
é, = —de, + re, + fe, 
the Lyapunov function and its derivative are yields Eqs. (15) and (16), respectively 
1 1 
V(e)=5 ener = 5 (er + ef + ef +e?) (15) 


V(e) = —ae? — qe? — ce? — de? < 0 (16) 


since V(e) is a positive function and V(e) is negative. So, the response of system (7) is PS with the drive 
system (6) asymptotically. The proof is complete. 
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4.3. The Controllers at scaling factor Y s; = —1 
For all scaling factor are taken the values -1, the error of PS between the system (6) and the system 
(7) is given by: 


é, = a(ez — e1) — fey + uy 


C2 = —Gez + 1X3 + €3X1 + (Y1 — X1) (V3 — x3) + u2 (17) 
È3 = —Ce3 — 61X2 — 2X — (Y1 — X1) (V2 — x2) + 2b + uz 
é4 = —de, + rez + u, 


Theorem 3. Choose the controller U; as: 


Uy = fe, — €2X3 + €3X2 


Uz = —ae; — re, — (x1 — y1) (x3 — Y3) (18) 
uz = W1 — %1)(Y2 — x2) — 2b 
U4 = 0 


The error dynamic system (17) can be controlled. 
Proof: With this control, (17) can be rewritten as: 


€, = a(ez — e1) — e2X3 + €3X2 

È> = —Gez — re, — ae + €1X3 + €3X1 
È3 = —C€3 — 1X2 — 2X1 

e4 = —de, + rez 


(19) 


The time derivative of the Lyapunov function is: 
V(e) = —ae? — qe? — ce? — dez < 0 


which is negative definite So, V(e) < 0. Therefore, PS of the two systems can be achieved simultaneously. 


5. HYBRID PROJECTIVE SYNCHRONIZATION (HPS) 
If at least one of scaling factor is different, this phenomenon is called HPS. Herein, two cases are 
considered as 
- gs, =1, 8. = 2, S3 = 3, s, =4 
= Seg hs 








5.1. The Controllers at scaling factor S4 = 1,52 = 2,53 = 3,5, = 4 
If the matrix S is chosen as S = diag(1, 2, 3, 4), i.e. 


y] pm 0 0 Qm 

_ |2 0 2 0 oļj* 

e= Ty3} 7 fa a (20) 
SERE 


Va X4 


won 


According to the (20), the error of HPS system between the system (6) and the system (7), is given as 


é = a(e, — e1) — fe, + ax, — 3fx,t uy, 


ê2 = —Gez + €1¥3 + (¥3—2%3)x1 + Ug (21) 
êz = —Ce3 — e2y1 — (2y1 — 3x1)xX2 — 2b + uz 
é4 = —de, + rez — 2rxX, + u, 


Theorem 4. If the following controller is designed as: 


Uy = —AxX2 + 3fxX4 — e2V3 
1 
uz = —5qe, + (2x3 — y3) x1 + €3¥1 — 500 Des 


uz = 2b + (2y1 = 3x4) x2 
u4 = 206re, + 2rx 


(22) 
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Then the system (21) will be controlled. 
Proof: Insert above control in (21), we get: 


ê, = a(ez — e1) — f e4 — €2Y3 
; 1 
ê, = —qez — 5qe, — 500 bea + e1Y3 + C34 


= —CEz — C21 
é, = —de, + re, + 206re, 


(23) 


D- 
w 
| 


The time derivative of the Lyapunov function is: 
V(e) = —ae? — qe? — ce? — de? + (a — 5q)e,e2 + (r — =b) eze, + (206r — f)e,e, (24) 


So, V(e) is negative definite, the system (21) was controlled based on control system (22). 





5.2. The Controllers at scaling factor s4 = 1,s2 = —1, S3 = 1,s4 = —1 
If the matrix S is chosen as S = diag(1, —1, 1, —1), i.e. 


Yı 1 0 0 07/% 
y2 0 -1 o0 0|% 
' H f 0 0 

y. 0 0 


S 


e= 


= 


0 —111%,4 
According to the above equation, the error of hybrid projective synchronization system between the system 
(6) and the system (7), is given as 


é, = a(en — e1) — fe, — ax, + fx, + uy 
êz = —qez + e3y1 + (V1 tx) x3 + U2 


; 25 
Èz = —Ce3 — e2y1 + (yı + X1)X2 + Ug Ca 
é, = —de, + re, + u, 
Theorem 5. If design the following controller (26): 
Uy = aX, — f X, 
Uz = —(x1 + y1)x3 — ae, — re, (26) 
uz = (—X, — y1)x2 
u, = fe 
Then the system (25) will be controlled. 
Proof: The time derivative of the Lyapunov function is: 
V(e) = —ae? — qe? — ce? — de? (27) 


So, V(e) is negative definite, the system (25) was controlled based on control (26). 


6. NUMERICAL SIMULATION 
For simulation, the MATLAB version R2017a is adopted to solve the differential equation 
of controlled error dynamical system (8), system (12) and system (17) for PS and controlled error dynamical 
system (21), system (25) for HPS based on fourth-order Runge-Kutta scheme with time step h = 0.01 
and the and the initial values of the drive system and the response system are following (3.2, 8.5, 3.5, 2.0) 
and (—3. 2, —8. 5, —3. 5, —2.0) respectively. We choose the parameters a = 5,b = 20,c = 1,d = 0.1,f = 
20.6, r = 0.1 andq = 1. 
— For scaling factor V s; = 3. Figure 2 show the PS of the systems (6) and (7) with control (9). 
— For scaling factor V s; = 1. Figure 3 show the PS of the systems (6) and (7) with control (13). 
— For scaling factor V s; = —1. Figure 4 show the PS of the systems (6) and (4) with control (18). 
— For scaling factor s4 = 1, s3 = 2, S3 = 3, s, = 4. Figure 5 show the HPS of the systems (1) and (4) 
with control (22). 
— For scaling factor s;3 =—1, S24 = 1. Figure 6 show the HPS of the systems (1) and (4) with 
control (26). 
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Figure 2. The PS for the state variables with control (9) at scaling factors V s; = 3 






























































Time(sec) Time(sec) 


Figure 3. The PS for the state variables with control (13) at scaling factors V s; = 1 
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Figure 4. The PS for the state variables with control (18) at scaling factors Y s; = —1 
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Figure 5. The HPS for the state variables with control (22) at scaling factors s; = 1, sz = 2, S3 = 3,s4 = 4 
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Figure 6. The HPS for the state variables with control (26) at scaling factors s; = +1 


7. CONCLUSION 

Based on the scaling matrix S, two types of synchronization phenomena are achieved, namely PS 
and HPS. Three error systems of PS and two error systems of HPS with controller have been proposed 
for obtaining PS and HPS between two identical 4-D hyperchaotic systems with unknown parameters based 
on Lyapunov's method and the nonlinear control strategy. Certainly, the projective synchronization, 
were achieved CS, AS as well as PS via this phenomenon. Whereas, the HPS, was achieved HS. 
The effectiveness of these proposed control strategies was validated by numerical simulation results. 
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